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BIT-LOADING IN MULTICARRIER COMMUNICATION SYSTEMS IN 
THE PRESENCE OF AN ASYMMETRIC, CORRELATED GAUSSIAN 
NOISE SOURCES 

TECHNICAL FIELD 

[0001] This disclosure relates generally to communication systems, and more 
particularly to multicarrier communication in the presence of noise and other 
sources of error. 

BACKGROUND 

[0002] A Discrete Multi-Tone (DMT) communication system carries 
information from a transmitter to a receiver over a number of tones. The tones 
are also commonly referred to as sub-carriers or sub-channels. The transmitter 
uses a modulation method in which the available bandwidth of a communication 
channel, such as twisted-pair copper media, is divided into these numerous sub 
channels. 

[0003] In the receiver, the data for each sub channel is typically extracted 
from the time-domain data by taking the Fourier transform of a block of samples 
from the multi-tone signal. There are various sources of interference and noise 
in a DMT system that may corrupt the information signal on each tone as it 
travels through the communication channel and is decoded at the receiver. 
Because of this signal corruption, the transmitted data may be retrieved 
erroneously by the receiver. In order to ensure a reliable communication 
between transmitter and receiver, each tone may carry a limited number of data 
bits. The number of data bits or the amount of information that a tone carries 
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may vary from tone to tone and depends on the relative power of the information 
and the corrupting signals on that particular tone. 

[0004] A reliable communication system is typically defined as a system in 
which the probability of an erroneously detected data bit by the receiver is 
always less than a target value. The aggregate sources of corruption associated 
with each tone are commonly modeled as a single additive noise source with 
symmetric Gaussian distribution that is added to the information signal on that 
tone. Under these assumptions, the signal-to-noise power ratio (SNR) becomes 
a significant factor in determining the maximum number of data bits a tone can 
carry reliably within a target bit rate error. 

[0005] The direct relationship between the SNR and the bit rate is based on 
the key assumption of the symmetry of the Gaussian noise source. However, 
this assumption may not be completely valid in many practical situations. There 
are certain types of noise disturbers that cause asymmetry in the distribution of 
the real and imaginary components of the noise. These noise sources can, for 
instance, be due to the inter-channel interference from the Nyquist tone - a sub- 
channel that only carry real component of the signal, or they can be due to ISI of 
background noise on stop-band frequencies. The effect of these sources 
become more complicated in a multicarrier system because the Frequency- 
Domain Equalizer rotates the signal and correlates the real and imaginary 
component of noise. The Frequency-Domain Equalizer also assumes that a 
symmetric Gaussian noise source is present in the background noise. With such 
a complex source of impairment, SNR alone cannot determine the reliable bit 
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rate. In the presence of an asymmetric Gaussian noise source, a typical bit- 
loading based on the assumption of symmetry may lead to an actual higher error 
rate than the target bit rate error. 
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SUMMARY 

[0006] In general, various methods and apparatus are described that use a 
multi-tone receiver. The multi-tone receiver detects data in a multiple tone 
signal. The receiver has a detector module to measure a noise power level 
present in the system and to detect for an asymmetric Gaussian noise source in 
the background noise. An equivalent noise power is obtained by applying a 
compensating gain factor to the asymmetrical noise source. The gain factor is 
used when the detector indicates that the asymmetric Gaussian noise source 
exists in the background noise. Any bit-loading algorithm that is based on a 
symmetric Guassian noise source assumption can also be used with 
asymmetrical Guassian noise sources if this gain factor is applied. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0007] Figure 1 illustrates a block diagram of an embodiment of a Discrete 
Multi-Tone (DMT) communication system containing a multi-tone receiver using 
a detector to detect whether an asymmetric Gaussian noise source is present in 
the background noise of a channel. 

Figures 2a-2b illustrate examples of scatter plots of noise error over 
time for a symmetric noise source and an asymmetric noise source. 

Figure 2c illustrates an approximate outlined representation of the 
noise distribution of the asymmetrical Gaussian noise source over time 
illustrated in figure 2b. 

Figure 3 illustrates a model for determining the equivalent total noise 
power of the asymmetrical, correlated Gaussian noise source. 

Figures 4a and 4b illustrates an embodiment a method of determining 
a signal to noise ratio in the presence of an asymmetric Gaussian correlated 
noise source. 
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DETAILED DESCRIPTION 

[0008] In the following detailed description of embodiments of the invention, 
reference is made to the accompanying drawings in which like references 
indicate similar elements, and in which, by way of illustration, specific 
embodiments in which the invention may be practiced are shown. These 
embodiments are understood that logical, mechanical, electrical, functional and 
other changes may be made without departing from the scope of the present 
invention. The following detailed description is, therefore, not to be taken in a 
limiting sense, and the scope of the present invention is defined only by the 
appended claims. 

[0009] In general, various methods and apparatus are described that use a 
multi-tone receiver, such as a DSL modem. The multi-tone receiver detects data 
in a multiple tone signal. The receiver has a detector module to measure a 
noise power level present in the system and to detect for an asymmetric 
Gaussian noise source in the background noise. An equivalent noise power is 
obtained by applying a compensating gain factor to the asymmetrical noise 
source. The gain factor is used when the detector indicates that the asymmetric 
Gaussian noise source exists in the background noise. Any bit-loading algorithm 
that is based on a symmetric Guassian noise source assumption can also be 
used with asymmetrical Guassian noise sources if this gain factor is applied. 
[0010] The equivalent noise power algorithm may alter the measured noise 
power level to model a noise power level of a symmetric noise source with a 
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radius equal to a major axis of the shape of the distribution of an asymmetric 
Gaussian noise error in a scatter plot. 

[0011] Figure 1 illustrates a block diagram of an embodiment of a Discrete 
Multi-Tone (DMT) communication system containing a multi-tone receiver using 
a detector to detect whether an asymmetric Gaussian noise source is present in 
the background noise of a channel. The DMT communication system carries 
information from a transmitter 102 through a channel .107 (such as a telephone 
line) to a receiver 104, such as a DSL modem, over a number of sub-carriers i.e. 
tones. In an Asymmetric Digital Subscriber Line (ADSL) system, each tone may 
be modulated using a Quadrature Amplitude Modulation (QAM) scheme. The 
transmitter 102 at a central office transmits the multicarrier signal (e.g. a DMT 
signal) to the receiver 104 using the QAM scheme. As the multicarrier signal 
travels to the receiver 104 through a transmission channel 107, the signal is 
potentially corrupted by various noise sources 106. 

[0012] The receiver 104 may contain various modules such as a Fast Fourier 
Transform block 110, block of Filters 1 12, a Total Noise Power Measurement 
block 114, Signal Power Measurement block 1 16, a Gain Factor block 1 18, as 
well as other modules. 

[0013] In the receiver 104, the data for each tone/sub-channel is typically 
extracted from the time-domain data by taking the Fourier transform of a block of 
samples from the multi-tone signal. The Fast Fourier Transform block 110 
receives the output of a block of filters 112. The Fast Fourier Transform block 
1 1 0 transforms the data samples of the multi-tone signal from the time-domain 
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to the frequency-domain, such that a stream of data for each sub-carrier may be 
output from the Fast Fourier Transform block 110. Essentially, the Fast Fourier 
Transform block 110 acts as a demodulator to separate data corresponding to 
each tone in the multiple tone signal. In one embodiment, processing of each 
sub-carrier may be performed in parallel or in series. The Fast Fourier 
Transform (FFT) block 110 may sample a sine and cosine of the amplitude of a 
tone over time to create the time domain data. The Fourier transform correlates 
the time domain data of the tone to the actual sine and cosine of the amplitude 
of the tone over time. 

[0014] For each particular sub-carrier of the multicarrier signal, the Total 
Noise Power Measurement block 1 14 measures the power level of total noise for 
that sub-carrier. 

[0015] The Total Noise Power Measurement block 1 14 measures noise 
present in the system by comparing the mean difference between the values of 
the received data to a finite set of expected data points that potentially could be 
received. The noise in the signal may be detected by determining the distance 
between the amplitude of the transmitted tone (at a given frequency and 
amplitude level) and the amplitude of the sine term and cosine term of the 
received tone to determine the magnitude of the error signal for that tone at that 
time. The noise present causes the error between the expected known value 
and the actual received value. The Total Noise Power Measurement block 114 
may contain a detector 115. 
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[0016] The detector 115 detects whether asymmetric noise is present in the 
background noise over time. The detector 115 module generates a scatter plot 
of noise error over time. The detector 115 analyses a shape of the distribution 
of the noise error in the scatter plot. The detector 115 may measure the x axis 
and y-axis of the shape of the distribution of the noise error. If an asymmetric 
noise source is present, then the detector 1 15 triggers the equivalent noise 
power algorithm to alter the measured noise power level to model a noise power 
level of a symmetric noise source with a radius equal to the major axis of the 
shape of the distribution of the asymmetric Gaussian noise error in the scatter 
plot. 

[0017] The error between the receive and transmit signal is usually depicted 
in a 2D scatter plot of error samples for each tone. The scatter plot has a 
horizontal x-axis and a vertical y-axis. The scatter plot illustrates the difference 
in amplitude between the known value of the transmitted data signal and the 
actual received signal over time. The Total Noise Power Measurement block 
114 may store a finite set of expected data points that potentially could be 
received and compares the actual received amplitude value of the tone to these 
expected amplitude values. In the scatter plot, each point corresponds to a data 
point and represents the error of the in-phase and perpendicular components of 
the carrier. 

[0018] Figures 2a-2b illustrate examples of such scatter plots of noise error 
over time for a symmetric noise source and an asymmetric noise source. Plot 
202 illustrates a scatter plot of the error signal when the sole source of error is a 
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symmetric Gaussian noise source. The shape of the noise distribution over time 
is roughly a circle. Thus, the diameter of the symmetric Gaussian noise source 
in the horizontal x-axis is approximately the same as in the vertical y-axis. In 
contrast, the scatter plot for the error from an asymmetric and correlated noise 
source has a more complex characteristic, as illustrated in plot 204. The 
asymmetry has caused the elongation of distribution in one direction and the 
correlation has caused the rotation of primary axis. The asymmetric shape of 
the error has an example elliptical shape having a major axis and a minor axis. 
The major axis of the elliptical noise distribution of the asymmetrical Gaussian 
noise source over time has a significantly greater diameter than the diameter 
along the minor axis. 

[0019] Figure 2c illustrates an approximate outlined representation of the 
noise distribution of the asymmetrical Gaussian noise source over time 
illustrated in figure 2b. Overall, the shape of the noise distribution of the 
asymmetrical Gaussian noise source 206 is elliptical and the major axis 208 may 
have a significantly greater diameter than the minor axis 210. 
[0020] Therefore, a simple symmetric Gaussian noise distribution may not 
accurately model the effects of the background noise. Specifically, using a 
simple power measurement for Signal to Noise Ratio (SNR) calculation 
underestimates the effect of noise and results in a higher bit-error rate. 
However, a more accurate model for background noise may take into account 
the long diameter of the asymmetric noise as the diameter of an equivalent 
symmetric noise source. 
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[0021] Referring to figure 1 , the detector 1 1 5 may trigger the Total Noise 
Power Measurement module 1 14 to use an equivalent noise power algorithm to 
compensate for the power level of an asymmetric Gaussian noise source in the 
background noise. The output of the Total Noise Power Measurement module 
1 1 4 feeds the gain module 118. 

[0022] The gain module 1 1 8 may determine total noise power level for each 
tone in the multi-tone signal based upon the equivalent noise power algorithm. 
The gain module uses the equivalent noise power algorithm to compensate the 
noise power level if the detector 115 indicates that an asymmetric Gaussian 
noise source exists in the background noise. Gain Factor block 118 also 
determines a Gain Factor associated with the asymmetric noise source. The 
equivalent noise power is input into a Signal-to-Noise Ratio (SNR) block 122. 
The Signal Power Measurement block 1 16 measures the signal power for the 
sub-carrier, and inputs the result into the SNR block 122. The SNR block 
determines a signal-to-noise ratio, which is used to determine bit loading 124 for 
all sub-carriers. 

[0023] As discussed, the detector 1 1 5 may rely on the fact that the power of 
an asymmetric noise source varies significantly along different axis. The power 
level along the strongest and weakest axis can be calculated from the 
parameters of the noise model. 

[0024] The general form of the probability distribution function of a two- 
dimensional, zero-mean Gaussian noise source is expressed as: 
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P(x, y) = 



1 



2nJ 




[0025] where o x and a y are the standard deviation along the two main axes 

and is the correlation between the components of noise along the two axis. 

The general Gaussian noise source is completely defined with these 3 
parameters. Note, a is the root mean square value and thus the standard 
deviation. 

[0026] The total power level of noise may be calculated as: 



for an uncorrelated noise source is zero. Therefore, a symmetric and 

uncorrected noise source can be defined by only a single parameter. In such 
case, the power along all axes is identical for a symmetric noise source. 
Therefore, for a symmetric and uncorrelated Guassian noise source: 



p = 
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g + a 

x y 



(2) 



[0027] For a symmetric Gaussian noise source, a x and a y are identical and 



a 



X 



= a 



y 



= G 



a 



= 0 



(3) 



P 



= 2a 



2 



Attorney Docket No.: 6491.P059 



-12- 



[0028] For an asymmetric, correlated noise source, power is not identical 
along all axes. The power value along the strongest and weakest axes can be 
calculated from the parameters of noise model as: 

jr (g, 2 +<W(^ 2 -<) 2+4 < 

2 

— 2 

(4) 

[0029] where P and P are the power values along the strongest and 
weakest axis, respectively. In the example of elliptical noise error scatter plot of 
figure 2b, the major axis of the ellipse would be the strongest axis. 
[0030] From equation (4), if the power along the strongest and weakest axes • 
is calculated to be identical, then the noise source is symmetrical and 
uncorrelated. If the power along the strongest and weakest axes is calculated to 
be different, then the noise source is asymmetrical and possibly correlated. 
[0031] Given a set of noise measurement samples over time, one can 
estimate the parameters of the noise model in a maximum likelihood fashion as 
follows: 

° I = ~ 
a ) = ~ 

° xy = x y 



1 N 

= — y y 2 

1 N 

1 y 1 = 1 

(5) 
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[0032] where x, and y, are the two components of the f h noise measurement 
sample. The total number of measurement samples over time is assumed to be 
N. 

[0033] Equivalent Noise Power 

[0034] A bit-error occurs when the noise amplitude is large enough so that a 
received QAM constellation point crosses the so-called decision boundary and is 
decoded as a neighboring point. The decision boundaries normally cross 
midway between the adjacent constellation points. For a symmetric Gaussian 
noise source, the ratio of the minimum-distance of constellation points to the 
power level of noise determines the bit-error rate. For a fixed power level of 
noise, the higher the minimum distance is, the lower the probability of error is. 
Bit loading determines the number of bits a given tone may carry. When 
maintaining a target error rate for a given power level of noise, the constellation 
size may be chosen such that its minimum-distance is 

d q = oc P (6 ) 

[0035] where d G denotes the minimum-distance when the noise source is 

purely Gaussian, P is the power of the symmetric Gaussian noise source and 
a is a constant factor that depends only on the error rate and the channel 
coding scheme. 

[0036] When the noise is not symmetric the error rate may be primarily 
dominated by the power of noise in its strongest axis. The outlying points on the 
strongest axis cause the most amount of error and are located the farthest from 
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the center of the constellation as depicted by the scatter plot. Therefore, for the 
purpose of maintaining the probability of error at some target value, an 
equivalent noise power estimation can assume conservatively that the 
asymmetrical noise source is equivalent to a symmetric noise source with a 
power per axes (a 2 ) equivalent to the power of the original noise source along 
the strongest axis. With this assumption, the power of the equivalent symmetric 
Gaussian noise source can be written as: 



[0037] This corresponds to a boost in measured noise power by the following ■ 



[0038] Therefore, for the purpose of bit-error rate analysis and bit-loading, an 
asymmetric, correlated noise source can be treated as a symmetric Gaussian 
noise source with the additional amplified power applied above. 
[0039] Figure 3 illustrates a model for determining the equivalent total noise 
power of the asymmetrical, correlated Gaussian noise source. The power of the 
noise is determined in the strongest axis and weakest axis. The summed 
powered is then multiplied by a gain factor to determine equivalent noise power 



(7) 



factor: 




1 + 




(8) 
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level of that noise when an asymmetric noise source is present in the 
background noise. Note that the gain factor, G a , is calculated to be unity for a 

symmetric, uncorrelated noise source. The equivalent total noise power from the 
gain module is input into a Signal-to-Noise Ratio calculation. 
[0040] Figures 4a and 4b illustrates an embodiment a method of determining 
a signal to noise ratio in the presence of an asymmetric Gaussian correlated 
noise source. 

[0041] The method 400 illustrated in Figures 4a and 4b constitutes a 
computer program made up of computer-executable instructions illustrated as 
blocks (operations) from 402 until 414. Describing the method by reference to a 
flpw chart enables one skilled in the art to develop such programs including such 
instructions to carry out the methods on suitably configured computers (the 
processor of the computer executing the instructions from computer-readable 
media, including memory). The computer-executable instructions may be 
written in a computer programming language or may be embodied in firmware 
logic. If written in a programming language conforming to a recognized 
standard, such instructions can be executed on a variety of hardware platforms 
and for interface to a variety of operating systems. In addition, the computer 
program is not described with reference to any particular programming 
language. It will be appreciated that a variety of programming languages may 
be used to implement the teachings of the invention as described herein. 
Furthermore, it is common in the art to speak of software, in one form or another 
(e.g., program, procedure, process, application, module, logic...), as taking an 
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action or causing a result. Such expressions are merely a shorthand way of 
saying that execution of the software by a computer causes the processor of the 
computer to perform an action or produce a result. It will be appreciated that 
more or fewer processes may be incorporated into the method illustrated in 
Figures 4a and 4b without departing from the scope of the invention and that no 
particular order is implied by the arrangement of blocks shown and described 
herein. For example, in one embodiment, processes represented by blocks may 
be performed in parallel. 

[0042] Overall, the equivalent noise power algorithm may automatically 
calculate the equivalent noise power in a multi-tone communication system in 
the presence of an asymmetric, correlated Gaussian noise source. The error 
samples are measured over a finite time interval of N samples for each sub- 
channel, The two main components of the f h measurement are labeled x, and y h 
The measurement is updated regularly after a new set of error samples are 
collected. 

[0043] In block 402, a receiver may determine a total signal power present for 
that tone over time. 

[0044] In block 404, a receiver may measure a power level of noise for each 
tone in a multi-tone signal. The equivalent noise power algorithm may calculate 
the average total noise power with the measured power level of noise. The 
equivalent noise power algorithm may determine the x axis amplitude of an error 
signal of the tone over time. The equivalent noise power algorithm may determine 
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the y axis amplitude of the error signal of the tone over time. The average total noise 
power may be calculated as: 



[0045] As discussed, the error samples are measured over a finite time 
interval of N samples for each sub channel. The two main components of the f h 
measurement are labeled x, and y,. The measurement is updated regularly after 
a new set of error samples are collected. 

[0046] In block 406, the noise power measurement algorithm may determine 
a Gaussian noise power level in the signal and if the noise source is generating 
an asymmetric pattern of noise. If an asymmetric, correlated noise source is 
present, then the average value of the x axis amplitude is significantly different 
than the average value of the y axis amplitude. The average difference noise 
power may be calculated as: 



[0047] Note, if the diameter of the noise distribution plot in the strongest axis 
is the same as the diameter of the noise distribution plot in the weakest axis, 
then the average power difference will calculate to be approximately zero. This 
indicates that a symmetric noise source is present and allows the equivalent 
noise power algorithm to be applied to solely symmetric noise sources to 
achieve the same results as the general form of the Gaussian noise model (2). 





yf) 
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[0048] The equivalent noise power algorithm may calculate the average 
cross-correlation of how the strongest axis relates to the weakest axis as 



[0049] The equivalent noise power algorithm determines the average values 
for the x axis amplitude, the y axis amplitude, and the correlation between the x 
axis amplitude and the y axis amplitude over time. 

[0050] In block 408, the equivalent noise power algorithm may calculate a 
gain factor associated with the asymmetric noise pattern. The equivalent noise 
power algorithm may calculate the equivalent gain factor for noise present in the 
background when an asymmetric noise source is present as 



[0051] When the calculated gain factor G a is compared to the expected gain 

of (1), then if the gain factor is > 1 , then asymmetric noise is present in the 
background noise of the channel. 

[0052] In block 410, the equivalent noise power algorithm may apply the gain 
factor to the total noise power level to calculate an equivalent total noise power 
of the effective symmetric Gaussian noise present in the system for a particular 
sub-carrier. The equivalent noise power algorithm may calculate the equivalent 
noise power when an asymmetric noise source is present as 
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= G . P , 



[0053] The final value of equivalent total noise power, a) , can be used in 

any bit-loading algorithm that is designed specifically for Gaussian noise 
sources. 

[0054] In block 41 2, the receiver may determine a signal-to-noise ratio based 
on a signal power of the signal and the calculated equivalent total noise power. 
As discussed, the calculated equivalent total noise power may be supplied to 
determine the SNR of a sub carrier. Subsequently, the SNR based on the 
equivalent total noise power may be used to determine the bit loading for that 
sub-carrier with asymmetric noise present on that channel. The resulting bit . 
error rate remains at target value even if the noise source is asymmetrical and 
correlated. Thus, after the gain adjustment is made, the equivalent symmetric 
Gaussian noise power level may be supplied to a bit loading/error rate algorithm 
to be used in the algorithm when calculating the bit loading for that particular 
channel. 

[0055] In block 414, the method 500 may be implemented more robustly by 
introducing a hysteresis basis on the calculated Gain Factor value. The 
equivalent noise power algorithm may determine if the gain factor exceeds a 
threshold amount. If so, the equivalent noise power algorithm activates the 
asymmetric Gaussian noise compensation based on an initial threshold and 
deactivates the asymmetric Gaussian noise compensation based on a reset 
threshold. Note, the initial threshold is greater than the reset threshold. 
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[0056] The result is an improved user experience, since fluctuation between 
asymmetric noise compensation and non-compensation is minimized. By way of 
example, consider an embodiment of the equivalent noise power algorithm that 
has a state in which the Gain Factor has a value of unity, e.g. G a =1 ; in other 

words, asymmetric noise compensation is inactive. In one embodiment, a 
threshold value of about G fl =1 .1 is used to activate the asymmetric noise 

compensation for the signal. In other words, until a Gain Factor is calculated as 
having a value of about 1 .1 or greater, the Gain Factor used in calculating the 
equivalent noise power will remain at a value of unity, e.g. G a =1 . However, in 

one embodiment, once the Gain Factor is calculated to have a value of about 
G a =1 .1 or greater, the asymmetric noise compensation is activated, and the 

actual calculated Gain Factor value is used to calculate the equivalent noise 
power. In another embodiment, once in an asymmetric noise compensation 
state, a lower threshold is used for deactivation of asymmetric noise 
compensation. For example, in one embodiment, where asymmetric noise 
compensation is active (e.g. G fl >1), the compensation is deactivated (e.g. G a is 
assigned a value of unity) only when the determined Gain Factor has a value of 
about G a * 1 . 03. Thus, the algorithm may employ set and reset states when 
implementing the asymmetric noise compensation. 

[0057] Thus, in one embodiment, from a state in which the Gain Factor has a 
value of unity, a first threshold is used to control when asymmetric noise 
compensation is activated. Once activated, a second threshold that is lower 
than the first threshold is used to control when the asymmetric noise 
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compensation is deactivated. It will be appreciated that other values may be 
used as thresholds for hysteresis, and that other measurements may be used to 
control activation and deactivation of asymmetric noise compensation. 
[0058] Embodiments of the present invention are applicable types of DSL 
systems, such as, but not limited to, ADSL, Rate Adaptive DSL (RADSL), Very 
High Bit Rate DSL (VDSL or VHDSL), High Bit Rate DSL (HDSL), Symmetric 
DSL (SDSL), ISDN DSL (IDSL), and Orthogonal Frequency Division Multiplexing 
(OFDM), as well as communications systems using other modulation 
techniques. Embodiments of the present invention are applicable to 
communication systems employing carrier signals in general. 
[0059] The above description of illustrated embodiments of the invention, 
including what is described in the Abstract, is not intended to be exhaustive or to 
limit the invention to the precise forms or embodiments disclosed. While specific 
embodiments of, and examples for, the invention are described herein for 
illustrative purposes, various equivalent modifications are possible within the 
scope of the invention, as those skilled in the relevant art will recognize. These 
modifications can be made to embodiments of the invention in light of the above 
detailed description. The terms used in the following claims should not be 
construed to limit the invention to the specific embodiments disclosed in the 
specification and the claims. Rather, the scope of the invention is to be 
determined entirely by the following claims, which are to be construed in 
accordance with established doctrines of claim interpretation. 
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